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Abstract 

We analyze the lepton sector of a Left-Right Model based on the gauge group SU(2) l x SU(2)r x 
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symmetry, we keep the right-handed symmetry breaking scale relatively low, while simultaneously 
satisfying phenomenological constraints on the light neutrino masses. We take the right-handed 
scale to be of order 10's of TeV and perform a full numerical analysis of the model's parameter space, 
subject to experimental constraints on neutrino masses and mixings. The numerical procedure 
yields results for the right-handed neutrino masses and mixings and the various CP-violating 
phases. We also discuss phenomenological applications of the model to neutrinoless double beta 
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I. INTRODUCTION 



The Standard Model (SM) of particle physics has thus far provided an incredibly accurate 
description of all observed experimental data. Nevertheless, the SM is widely regarded as 
being a low-energy effective theory, with a limited range of applicability and predictability. 
New interactions must arrive at the energy scale of several TeV in order to explain such 
features of the SM as the quark and lepton mass hierarchy. 

An intriguing aspect of the SM is that its weak interaction sector represents the only 
known interaction that distinguishes between the right- and left-handed fermions. This 
aspect is addressed in a set of aesthetically-pleasing "left-right models" (LRMs) based on 
the gauge group SU(2)l_x SU(2)rX U{t), where the left- and right-handed fermion fields are 
treated symmetrically 

ft 3 13 BE IE 0- in such models, left-right symmetry is broken at 

some high scale, yielding a parity-violating Standard Model-like theory at low energies. In 
typical phenomenological studies of the LRM, the right-handed scale (i.e., the scale at which 
SU(2)r is broken) is assumed to be very high, of the order of 10 10 GeV. Such a high energy 
scale would render direct experimental verification of a LR-motivated scenario impossible in 
the near future. By way of contrast, more moderate values for the right-handed scale - in 
the range 20-50 TeV, say - could have observable consequences at experiments in the near 
future [E This is the energy range that we shall consider in this work. 

Recently, there has been a renewed interest in the LRM, in particular, due to the discovery 
of neutrino oscillations and to major advances in experimental studies of CP-violation in the 
quark sector. Of interest to us in this work is that the LRM provides a natural process for the 
suppression of neutrino masses through the seesaw mechanism. LRMs also offer additional 
sources of CP violation, coming both from the right-handed Cabibbo-Kobayashi-Maskawa 
(CKM) and Maki-Nakagawa-Sakata (MNS) matrices as well as from the Higgs sector of the 
theory EE| . The CP violation occuring in the leptonic sector of the model could in principle 
be of interest within the context of leptogenesis 11, EE EE EH EE EE E3|- ^ * s important 



to emphasize that the phases in the left-handed MNS matrix are currently unconstrained. 

As noted above, the LRM contains the natural possibility of implementing the seesaw 
mechanism for the generation of small neutrino masses. In the simplest version of the 
seesaw mechanism, only the SM-singlet right-handed neutrinos are initially allowed to obtain 
Majorana mass terms, resulting in the neutrino mass matrix 

where Mrr and M^r are Majorana and Dirac mass matrices, respectively, in flavour-space. 
This construction, with a block of "zeros" where the left-handed Majorana mass terms would 
go, leads to what is sometimes called the "Type I" seesaw mechanism. It can be generated 
in the LRM. An approximate block-diagonalization of Eq. ([TJ, assuming that the elements 
in Mrr are much larger than those in Mlr, leads to the standard seesaw expression for the 
light neutrino mass matrix, M v ~ —MlrM^ r MJ^ r = M£. 

To implement the seesaw mechanism in the LRM, one introduces a right-handed Higgs 
triplet field, Ar, into the theory. The field Ar serves a dual purpose - it breaks the left-right 
symmetry of the model at a high scale and it also couples to Majorana neutrino fields, giving 
rise to the right-handed Majorana mass matrix Mrr required for the seesaw mechanism. 
Mrr is proportional to the right-handed symmetry breaking scale, making it naturally 
large. Left-right symmetry also requires the existence of a left-handed Higgs triplet field, 
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Al. The neutral components of the left- and right-handed triplet fields both generically 
obtain a vacuum expectation value (VEV) under spontaneous symmetry breaking, (51 R ) = 

VL,Re t6h ' R / \/2. LRMs also typically contain a bidoublet field, 0, with VEVs at the weak 
scale. The role played by <fi is similar to that played by the usual Higgs doublet of the SM. 

In contrast with the situation in the Type I seesaw mechanism, LRMs generically also 
contain a non-zero left-handed Majorana mass matrix for neutrinos, Mll- This mass matrix 
is proportional to v L (the VEV associated with the left-handed Higgs triplet) and vl need 
not be zero. There is, in fact, a seesaw-like relation among the VEVs vl and v R in the LRM 
that arises due to terms in the Higgs potential that couple 0, A L and A R (Tr(0A R 0^A^), for 
example). If all dimensionless coefficients in the Higgs potential are of order unity, one finds 
that vl ~ k 2 /vR, where k is a dimensionful quantity of order the weak scale [ly, [3, fl9l | . 
The situation v l = (and hence Mll = 0) can be obtained by dropping the offending terms 
from the Higgs potential, although attempts to disallow the terms using a symmetry meet 
with difficulties [19]. Inclusion of a non-zero matrix M LL in the "zero" block of Eq. ((TJ (see 
Eq. (|15j) below) leads to the "Type II" seesaw mechanism [lil l2o| . yielding the following 
approximate mass matrix for light neutrinos, 

M v ~ M{ L - M LR M- R \M T LR ee M l J + Ml (2) 

with Mll generically of order vl ~ k 2 /v R and M RR ~ v R . Experimentally, the terms in M v 
must be at most of order about 0.1 eV. If the largest terms in Mlr are taken to be of order 
m T , and if Mll does not undergo any further suppression, then the first term in Eq. (J2J) 
dominates and sets the minimum scale v R that is phenomenologically viable. Taking k to 
be of order the weak scale we find that v R would need to be at least of order 10 14 GeV in 
this case, ruling out any possibility of observing LR-induced effects at collider experiments. 
The analogous lower bound coming from the second term alone (Type I seesaw) is of order 
10 10 GeV if one assumes M L r ~ m T . 

A few approaches have been suggested for reducing the right-handed scale v R while si- 
multaneously satisfying phenomenological constraints within the context of Type II models. 
One approach is to suppress Mll by separating the parity and gauge symmetry breaking 
scales. For instance, one can introduce a pseudoscalar Higgs field 77 that acquires a vacuum 
expectation value (thereby breaking parity) at a very high scale, while allowing the right- 
handed gauge symmetry to be broken at a much lower scale [2lJ, |22| This approach could lead 
to interesting phenomenology. In such a model, M L l ~ k 2 v R / (r]) 2 , so that the left-handed 
Majorana mass terms could be of an acceptable size provided that (rj) were sufficiently high. 

Another promising approach for bringing v R down to a potentially observable scale is to 
introduce an extra £7(1) horizontal symmetry that is broken by a small parameter e [ToL I23L 
0- HH . In this approach, each field in the model is assigned a charge under the horizontal 
symmetry. Yukawa couplings and the dimensionless coefficients in the Higgs potential are 
then suppressed by various powers of e, with the powers depending on combinations of 
charge assignments. In addition to providing a nice dynamical mechanism for producing the 
observed hierarchies in the charged lepton masses, this approach also allows for a significant 
reduction in the right-handed scale v R . As shown in Refs. (HI an appropriate choice 
of charge assignments leads to a suppression of vl, thereby suppressing M L l- There is a 
similar suppression of the Yukawa couplings involved in Mlr, the net effect of which is to 
loosen the stringent lower bounds on vr. 

In this work we perform an indepth numerical study of neutrinos in a left-right model 
that is supplemented by such a broken U(l) horizontal symmetry. The goal of the paper is 
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to determine the parameter space of the model already constrained by the neutrino mass 
and mixing measurements. We focus on the case in which the right-handed symmetry 
breaking scale is only "moderately" large (20-50 TeV). In the horizontal symmetry scheme, 
right-handed scales as low as 20 TeV can in fact lead to results consistent with neutrino 
phenomenology. The numerical procedure also yields results for the right-handed neutrino 
masses and mixings and the various CP-violating phases. Throughout this work we assume 
that there are three generations of light neutrinos and ignore the possibility of light sterile 
neutrinos. We also assume the "normal" ordering of neutrino masses, m Ul < m U2 < m Us . 

The outline of the remainder of the paper is as follows. In Sec. |H]we describe the LRM and 
establish our notation. Section IIHI contains our numerical results for the allowed parameter 
space of the model, as well as a discussion of some phenomenological applications such as 
lepton-flavor violating transitions and leptogenesis. We conclude with a brief discussion 
in Sec. II V I The Appendix contains approximate expressions for the right-handed neutrino 
masses and mixings as well as a particular case study. 



II. THE MODEL 




A. Mass and MNS matrices 

The LRM is based on the gauge group SU(2) L x SU(2) R x U(1)b-l- We will consider 
a minimal version of the model, whose Higgs sector contains a bidoublet Higgs boson field 
(p ~ (2, 2, 0) and two triplet Higgs boson fields A L ~ (3, 1, 2) and A R ~ (1, 3, 2). As noted 
in the Introduction, A R is used to break the gauge symmetry down to SU(2)l x U{1)y and 
(ft is used to break it down to £/(l) cm . The various Higgs boson fields may be parametrized 
as follows, 

The left- and right-handed lepton fields transform as doublets under SU(2)l and SU(2) R , 
respectively, and are given by 

Vk,H = ( f' R ) , (4) 

where i is a generation index and where the primes denote that the fields are gauge eigen- 
states. In addition to the gauge symmetry, it is common to impose an extra left-right parity 
symmetry 0], demanding invariance under 

^'iL^^'iw 4>^<t>\ A L ^A R . (5) 

The lepton Yukawa couplings that are consistent with the gauge and parity symmetries 
discussed above are [l^] 

- £ Yukawa = (Gij<P + Hy$) i>' jR + \Fij (^aCr 2 A L iP' jL + ^ R Cr 2 A R ^' jR ) + h.c. , (6) 
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where <ft — t~ 2 <J)*t 2 and where C = ry 2 7° is the charge conjugation matrix. In order to respect 
the parity symmetry in Eq. (jSJ), the 3x3 matrices G and H must be Hermitian. The matrix 
F is complex, in general, but may be taken to be symmetric without any loss in generality. 1 
Upon spontaneous symmetry breaking the Higgs boson fields acquire VEVs, which may 
be parametrized as 

with ki 7 2 an d v l,r rea l an d positive. Gauge rotations have been used to eliminate possible 
phases associated with k\ and vr Phenomenological constraints require that vr ^> 

ki, k 2 3> In this case ki and k 2 satisfy the constraint |27] 

k \ + k\ ~ ~ (246.2 GeV) 2 . (8) 

Also, it is natural to assume k 2 jk\ ~ mb/m t ; in the numerical work below we shall set 
k 2 /h = 3/181 as in Ref. @. 

The VEVs in Eq. (J7J) lead to Dirac mass terms for the neutrinos and charged leptons, 

- ^Dirac = ^ (GAu + ^fc 2 e" ia ) ^ + i= i£ (Gk 2 e ta + Hk x ) e' R + h.c. , (9) 

as well as Majorana mass terms for the neutrinos, 

" ^Majorana = ~^ F V ^ L v' L + I^FVrU'r) + h.C, (10) 



where u'£ R = Cv' LR . The mass matrix for the charged leptons is thus 

M e = ^= (Gk 2 e ia + Hh) , (11) 



V2 

which may be diagonalized by a biunitary transformation 

Mf ag = V?M t Vji , (12) 

where the elements in M^ iag are real and positive. 

Consideration of the neutrino mass matrix is slightly complicated by the fact that both 
Majorana and Dirac mass terms are present. Defining 

v " = {%) (13) 

allows the neutrino mass terms in the Lagrangian to be written as 

-£ m „ = ^W u MK c + h.c., (14) 



1 This follows because one may write, for example, vf L v^ L — v'^ L v[ L ^ where v'^ L = Cv' iL 
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where the Majorana and Dirac neutrino mass matrices have been incorporated into a single 
6x6 complex symmetric matrix 



with 

M, 



lr = 4= (Gh + Hk 2 e~ ia ) , M LL = l=Fv L e iB \ M RR = ^Fv R . (16) 

The 6x6 neutrino mass matrix Ai may be approximately block diagonalized into a light, 
mostly left-handed block, M u , and a heavy, mostly right-handed block, M R , 



m v o \ / 1 e\ T (m[ l m lr \ ( i ? 

M R ) ~ [-£ 1) \M T LR M RR ) \ -£ 1 

„ (M\ L - M LR M RR Ml R \ 
V M RR 



(17) 
(18) 



where £ = M RR Ml R . The corrections to Eqs. (|T7|) and (fTHj) are suppressed by £ and may 
be neglected if 1. (In our numerical work below, the magnitudes of the elements of 

£ are of order 1CT 5 or smaller, so the approximation is well-justified.) The above result for 
M v was quoted in Eq. (j2J) and is the general expression for the Type II seesaw mechanism. 
While the Dirac mass matrix M^ R is of the same order of magnitude as its counterpart for 
the charged leptons ("M/'), its contribution to M v is suppressed by £. 

The six physical neutrino states obtained by diagonalizing M. are all Majorana neutri- 
nos [26], a fact that is responsible for the symmetry of M. about the diagonal. The two 3x3 
blocks in Eq. (JTSj) are also symmetric and may each be diagonalized with a unitary matrix, 
yielding 

M dia g = y£ Mv yv* f ( 19 ) 
M dia g = V f MRRV » f (20) 

where the elements of the diagonal matrices are real and positive. 

The unitary matrices V[ R and V£ R used to diagonalize the charged and neutral lepton 
mass matrices may be combined to give the MNS matrices, 

QMNS = yiiyu (21) 

= V^Vr , (22) 

where the tildes indicate that the matrices may still be "rephased" to bring them into 
the conventional form. The rephasing procedure for the MNS matrices is accomplished by 
multiplying the expressions in Eqs. (|21|) and (|22|) on the left and right by diagonal phase 
matrices, 

U? N3 = B^ INS S L , (23) 
U^ NS = B^ INS S R , (24) 
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where 



B 



( e ipl 

e lp2 ! . S, 
\ e ip 



e im2n | , S R 




/ e im ^ 

e im * n | , (25) 
V e im ^ 
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with rrii and mf integers. The above rephasing procedure differs from its counterpart in the 
quark sector in two respects. In the first place, since the neutrinos are Majorana particles, the 
"phase" matrices Sl and Sr are actually only "sign" matrices, with factors of ±1 appearing 
along the diagonals. A second difference compared to the quark case is that Sl and Sr 
are distinct matrices, whereas in the quark case the analogous matrices are equal We 
may use the above results to write the charged current couplings in the Lagrangian in terms 
of the physical mass eigenstates. Ignoring terms further suppressed by £ in Eq. (JIBJ), we 
have [26] 

where 

"l,r = SIrV£ r v LjR , (27) 
e L , R = B^ R e' L>R . (28) 

It is in general possible to parametrize the rephased left-handed MNS matrix in terms of 
three non-removable CP-odd phases. One of these phases is analogous to the usual CKM 
phase in the left-handed CKM matrix. The other two phases are novel, compared to the 
quark sector, and their presence is due to the Majorana nature of the neutrinos. Attempts 
to remove these Majorana phases result in their appearing elsewhere in the theory (in the 
diagonalized neutrino masses, for example). A useful parameterization of the left-handed 
MNS matrix is ET 



u^ INS = u^\e 12 ,e 23 ,e 13 ,s L )A L , (29) 

where A L = diag(e* Ql / 2 , e ia2 ^ 2 , 1) and 
U {0 \0 12 , 9 23 , 9 13 , 6 L ) 

' c 12 c 13 s u c 13 s 13 e~ l&L \ 

= -S12C23 - c 12 s 23 s 13 e lSL c 12 c 23 - s 12 s 23 s 13 e lSL s 23 ci 3 . (30) 

\ s 12 s 23 — c l2C 23 Si 3 e l L —C\ 2 S 23 — Si2C23Si3e* L C23C13 J 

The phase 8l in the above expression is analogous to the usual CP-odd "Dirac" phase in 
the quark sector, while the phases ot\ and a 2 are Majorana phases. The right-handed MNS 
matrix contains six phases in general, three of which may be taken to be Majorana phases. 
A convenient parameterization is as follows, 

Ur S = 4u^\ef 2 ,eg,9f 3 ,8 R )AR, (31) 

where A R = diag(e ia ? /2 , e^ /2 , e ia * /2 ) and A e = diag(e iCl , e i?2 , 1). 

The left-handed MNS matrix has been probed through neutrino oscillation experiments, 
which have placed relatively tight constraints on the three mixing angles 9ij. The left-handed 
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phases 5l and a± t 2 have not as yet been constrained by experiment. Neutrinoless double 
beta decay experiments could well be used to probe combinations of the left-handed phases 
(depending on the ordering of the light neutrino masses - i.e., "normal" or "inverted" - and 
on the magnitudes of the masses). In fact, such experiments play a central role in neutrino 
physics, since the decays in question can only proceed if neutrinos are Majorana (as opposed 
to Dirac) particles. The amplitudes for such decays are proportional to mpp, the effective 
neutrino mass for neutrinoless double beta decay |26j], 



e Kr ) 



J m 3 



\m lC 2 12 c 2 13 + m 2 s 2 12 c 2 13 e iia2 - ai) + m 3 s 2 13 e- i(ai+2SL) 



(32) 



Future experiments could probe mpp at the O(10~ 2 eV) level (see, for example, Ref. j28j . 
as well as ). 2 Evidently could be a sensitive probe of the MNS phases if the mixing 
angles and masses were well known. 3 In the numerical work below we will calculate mpp for 
this model to determine prospects for future experiments. 



B. Simplification of the Yukawa Couplings 

It is often possible to simplify the Yukawa couplings in a model through unitary trans- 
formations that leave physical quantities, such as masses and mixings, unchanged. In the 
present case, unitary transformation on the Yukawa coupling matrices G, H and F may be 
used to reduce the number of parameters required to specify the model. As noted above, in 
order to satisfy the parity symmetry in Eq. (|5|). G and H must both be Hermitian. Further- 
more, F may be taken to be (complex) symmetric. For three generations of leptons, this 
means that 30 real parameters are required to specify the elements in the Yukawa matrices. 
In principle, several of these degrees of freedom are spurious and may be "rotated away" by 
an appropriate unitary rotation. To see this, note that the diagonalized mass matrices and 
MNS matrices are invariant under the rotations 

F -> X T FX, 

G -> X*GX, (33) 
H -»• X*HX, 

where X is a 3 x 3 unitary matrix. Furthermore, these rotations preserve the essential 
symmetries of the matrices, leaving G and H Hermitian and F symmetric. In principle, one 
could use X to diagonalize one of the three Yukawa matrices, significantly decreasing the 



2 There is controversial evidence of a non-zero neutrinoless double beta decay signal with mpf} of order 
0.5 eV 00- 

See also Refs. 

3 In our notation the expression for mpp contains the phase combinations ct2 — oc\ and ai+2#j,. It is possible 
to rephase the MNS matrix in such a way that the 1-3 element of Z^ -* is real and the two Majorana phases 
in Al occur in the 2-2 and 3-3 elements j^- I n that notation mpp depends only on the (two) Majorana 
phases contained in Al . We have verified that the relations between the phases used in the two approaches 
are such that one obtains the same physical value for mpp in either approach. 



number of parameters required to specify the model. Within the context of the horizontal 
symmetry scheme that we employ, however, the above transformations affect the scaling 
of the various terms. For this reason we do not diagonalize any of the Yukawa matrices, 
choosing instead to use a phase rotation in Eq. to remove one phase from F (arg(F 22 )) 
and two from H (arg(ifi 2 ) and arg(if 13 )). This reduces the number of parameters required to 
specify the Yukawa couplings to 27. In the numerical work below, a Monte Carlo algorithm 
is used to search the 2 7- dimensional parameter space to determine sets of parameters that 
are consistent with experimental constraints on the lepton masses and mixings. 



C. Model with a broken U(l) symmetry 

One attractive way to account for the observed hierarchies in the quark and lepton Yukawa 
couplings is to attribute them to a broken horizontal symmetry [i^, 24]|. In models with 



a broken horizontal symmetry, the various Yukawa couplings are suppressed by powers of 
one or more small parameters, where the powers are determined by the charges of the 
relevant fields under the horizontal symmetry group. Khasanov and Perez [25[ recently 
formulated a model that uses a broken horizontal U(l) symmetry to address two known 
problems that occur in the LRM if one attempts to take v r to be only moderately large 
(of order 20 TeV, say). The two problems are associated with the two terms appearing in 
Eq. (0) - as noted in the Introduction, both terms run into trouble for moderate values 
of vr unless they are suppressed in some manner. The first term in this expression is 
particularly troublesome - although it is somewhat suppressed due to the VEV seesaw, it 
is still far too large. For vr of order 20 TeV, minimization of the Higgs potential yields 
vl ~ k 2 /vR ~ (246.2 GeV) 2 /(20 TeV) ~ 0{l GeV) assuming all dimensionless coefficients 
in the Higgs potential to be of order unity (see Ref. [10j . for example). If the Yukawa matrix 
F in (JTfij) is of order unity, then Mll will be of order 1 GeV, approximately nine or ten 
orders of magnitude larger than the neutrino mass scale. The second term in Eq. (J2J) is 
also too large if vr is of order 20 TeV. Assuming F to be of order unity and the largest 
elements of M LR to be of order m T , we find that u —M LR M RR Ml R " generically has elements 
of order m 2 T /v R ~ (1.777 GeV) 2 /(20 TeV) ~ 0(0.1 MeV), which are still too large from a 
phenomenological point of view. 4 

A model with a broken horizontal U(l) symmetry offers a solution to both of the problems 
noted above. At high energies the model contains a new scalar S as well as several new 
heavy fermions. Most of the Yukawa terms in Eq. (jUJ) are not present in the high energy 
theory because they do not respect the U(l) symmetry. Instead, such terms descend from 
nonrenormalizable terms in the low energy effective theory obtained by integrating out the 
heavy fermions. As a result of this procedure, the Yukawa terms contain various powers of 
a small symmetry breaking parameter e = (S)/M, where M is the mass scale of the heavy 
fermions. The power of e for a given term in the Lagrangian is determined by the U(l) 



4 One could improve the situation by assuming that the Yukawa matrix G ~ 0. In that case, the largest 
elements in H are of order rrij/ki and we have Mlr ~ Hkz ~ m T x (fe/fci). As noted above, it is natural 
to assume fcg/fci ~ m^jmt [ij, m which case the largest elements in "— MlrM^Mlh" are of order 10's 
of eV for vr — 20 TeV. In the horizontal symmetry scheme that we employ below, G is in fact suppressed 
relative to H, leading to a similar result (see Eq. l|42|l and the discussion that follows). 
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charges of the fields coupled together in that term. The Yukawa couplings scale as follows, 



p.. _ p. \Q(A L )+Q(Li)+Q(Ll)\ 
1 ij — 1 ijt 



Gij = ^. e IQ(^)-Q(^)+W)l ; (34) 



where the quantities with the tildes are taken to be of order unity in magnitude. In our 



igmtucu 

numerical work we adopt the following charge assignments (see also Refs. EEQ) 

Q(A L ) = -Q(A R ) = -8, 

W) = -2, (35) 

Q(4' 2 < 3 ) = -Q(4 2 < 3 )= 6,4,3, 

yielding 
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(36) 



where coefficients of order unity have been omitted. For the purpose of our numerical work 
we set e = 0.3, as in Ref. a value that automatically gives charged lepton masses in the 
correct range. The Higgs potential of the low energy effective theory also contains terms that 
break the U(l) symmetry, leading to a suppression of many of the dimensionless coefficients 
in the Higgs potential. Reference contains a thorough discussion of the Higgs sector of 
the LRM with a broken horizontal symmetry. 5 With the charge assignments noted in (J33j) . 
minimization of the Higgs potential leads to the following expression for vl (l0|, 

v L = ie 20 k\/v R , (37) 

where 7 depends on various dimensionless coefficients in the Higgs potential and is generically 
of order unity. For 7 of order unity (and setting e = 0.3), we have 

Bi ~ •WPOTeV)' (38) 

which is phenomenologically viable. Thus the U(l) model successfully deals with the first 
of the two problems noted above. 

The U(l) model also deals successfully with the fact that the second term in Eq. (j2J) is 
generically too large. Before considering this term, let us examine the charged lepton mass 
matrix, Mg, given in Eq. (jllj) . To a good approximation, one may neglect the contribution 
of G to Mi, since this contribution is suppressed by a factor of approximately e 4 x (k 2 /ki) ~ 
O(10~ 4 ) relative to that of H . The situation is different for the neutrino Dirac mass matrix 



5 In that paper it was shown that a phenomenologically acceptable Higgs spectrum emerges if explicit CP vi- 
olation is allowed in the Higgs potential. This is to be contrasted with the case in which the Higgs potential 
is CP-invariant. In that case, non-negligible CP violation in the vacuum state is generically accompanied 
by non-SM-like neutral Higgs bosons at the weak scale with flavour non-diagonal couplings |35| . 
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Mlr, since the roles of k\ and k 2 are essentially reversed in this case. In fact, G and H 
contribute comparable amounts to Mlr (since fe/^i ~ C( e4 )) an d we find that Mlr = 
(Gh + Hk 2 e-' ia ) /y/2 ~ Afc x e 4 as an order of magnitude estimate. Noting that 



(39) 




(where a "1" denotes an element of order unity), we have 



ML 



M LR M RR M T LR 



k\ 



V2v R 
V2v R 



( 



,24 ^.22 ,2\ 
z t t 

,22 ^20 ,-21 



+ 



k\ 



,21 J2.1 



,22 



£ 24 £ 22 e 21 
,22 ^20 ^19 



,21 ^19 ^18 



,24 ,22 ,21 
z t t 

,22 c 20 ,19 
zee 

,21 A9 A8 



v R /(20 TeV) 



0.0006 0.007 0.02 
0.007 0.07 0.2 
0.02 0.2 0.8 



eV, 



(40) 



(41) 



(42) 



where the numerical values in the last line should be understood as being very approximate. 6 
The two terms in the above expression for M u contribute at approximately the same 
level and combine to yield neutrino masses that are of the correct order of magnitude. 
It is interesting to see how the horizontal symmetry model deals with the fact that the 
largest elements in the Type I seesaw part of M v (the second term) are generically of order 
w? T jvR ~ 0.1 MeV. The main suppression of such elements in the £7(1) model follows from 
the fact that the largest terms in Mlr are now of order e 4 m r , instead of m T , as noted in the 
discussion above Eq. ([3~9*]). A further suppression is due to the particular structures of F -1 
and M LR . 7 

The LRM with a broken £7(1) symmetry is thus able to reproduce the gross features of 
the lepton mass spectra, yielding the correct orders of magnitude for the charged lepton 
masses as well as an appropriate mass scale for the light neutrinos. In the following, we 
consider whether the model is able to accommodate the experimental values for the light 
neutrino mass-squared differences and mixing angles. In fact, there is potentially a difficulty 
in this regard, as was pointed out by Khasanov and Perez [25j - the 1-2 element in Eq. ()41|) is 
generically suppressed relative to the 2-2 element, indicating a possible difficulty in obtaining 
a large 1-2 mixing angle. Nevertheless, we shall show that it is in fact possible to satisfy 
the experimental constraints on all three mixing angles and on the mass-squared differences 
in this model. We offer some further comments on this issue in Appendix As is noted 
there, the numerical procedure favours neutrino mass matrices that have a quasi-degenerate 
2-3 block, with some or all elements in the block suppressed relative to Eq. 



6 Given the large exponents in this expression, one might worry about the sensitivity of these results to 
small deviations in the parameter e. While it is true that a small change in e would produce a larger effect 
in M„, the effect on M v could be partly compensated by adjusting vr. 

7 For example, one contribution to the 3-3 element in M v comes from the 3-3 elements of Mlr. and F~ l . 
While Mlr.33 is generically the largest element of Mlr, F^ 1 ~ e 2 , so the combined contribution is of 
order e 2 x (e 4 m T ) 2 /wfl. 
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III. NUMERICAL STUDY 

In this section we perform a numerical analysis of the model. The goal of the numerical 
work is to find sets of values for the various Higgs VEVs and Yukawa couplings such that 
the experimental constraints on the lepton masses and mixings are satisfied. Equation (J34j) 
expresses the three Yukawa matrices F, G and H in terms of rescaled (order unity) Yukawa 
couplings (Fij, etc.) multiplied by appropriate powers of e. Many of the Yukawa couplings 
are complex. Recalling that F is complex-symmetric and that G and H are both Hermitian, 
we define phases as follows, 



F 


= F = 


F ■ 


i8 F 

e ^ 


U>i), 




Gij 


~ G ii - 


Gij 


i8 G 

e ^ 


> 0, 


(43) 




= H* 3l = 


Hij 


ie H 

e « 


U > *)■ 





The diagonal elements of G and H are real, but possibly negative. As described in Sec. Ill 131 
unitary rotations may be used to simplify the Yukawa matrices without affecting the lepton 
masses or the MNS matrices. We use such rotations to eliminate one phase in F and two in 
H, setting 0$ 2 = 6^=6^ = 0. Thus, there are a total of 27 parameters used to describe the 
Yukawa matrices, nine of which are phases. In our numerical work, we allow the magnitudes 
of the scaled Yukawa couplings to be in the range zero to three and the phases to be in the 
range zero to 2ir. 

For the Higgs VEVs, we use Eq. (jHJ) to fix the sum k\ + k\ and take the ratio k 2 /ki 
to be 3/181 (as in Refs. 0,[[(3)- We consider two cases for vr, taking vr = 20 TeV and 
vr = 50 TeV. vl is defined through Eq. (J37|) . where we take 7 to be chosen randomly in 
the range zero to two. It remains to consider the phases of the Higgs VEVs, a and 9l (see 
Eq. (J7|)). Correlations between a, 9l and vl were studied in Ref. [lOj. Since the observed 
correlations were not very strong, we simply allow a and 6l to take any values in the range 
zero to 2n. Adding a, 9l and vl to the 27 Yukawa coupling parameters, we find that we have 
a total of 30 parameters to fix. This number exceeds the number of experimental constraints 
on the model, which come from the charged lepton masses (3), the neutrino mass-squared 
differences (2) and the neutrino mixing angles (3). 8 Clearly it will not be possible to fix 
the 30 "input parameters" uniquely. Nevertheless, a Monte Carlo approach can be used 
to find sets of input parameters that yield masses and mixings consistent with experiment. 
Once these experimental constraints have been satisfied, other quantities in the model - 
such as left-handed neutrino phases and right-handed neutrino masses and mixings - can be 
calculated. 

A. Monte Carlo Algorithm 

The Monte Carlo approach that we use is similar to that described in Ref. 0. A rough 
summary of the procedure is as follows. Sets of input parameters are chosen randomly 
and then used to form the various mass matrices. Diagonalization of these mass matrices 
yields theoretical values for the masses and mixing angles, which are then compared with 



We do not include the LSND results in our analysis. 
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TABLE I: Experimental constraints for neutrino masses and mixings used in the numerical work. 
The values adopted for y^ xp and <7j correspond to the 3a ranges given in Table 1 of Ref. |37|. We 
estimate the central values y^ xp by bisecting the 3a ranges. Normal mass ordering is assumed in 
the numerical work, so both of the mass-squared differences are taken to be positive. 

Quantity y exp ± ai 

Am^ = m \-m\ (8.15 ±0.95) x 1(T 5 eV 2 

Am§! =m\-m\ (2.35 ± 0.95) x 10~ 3 eV 2 

sin 2 6i2 0.305 ± 0.075 

sin 2 23 0.51 ±0.17 

sin 2 0i3 0.0235 ± 0.0235 



their experimental counterparts. Specifically, we calculate the charged lepton masses, the 
neutrino mass-squared differences and the squares of the sines of the mixing angles and 
compare these to the experimental values described in Table |U 9 A quantitative measure of 
the "goodness of fit" is provided by the quantity x 2 , 

8 |V X P _ v .)2 

where the sum runs over the five experimental constraints y^ xp ± a, in Table |1J as well as 
three constraints coming from the charged lepton masses. The associated values obtained 
numerically are denoted i/i. The Monte Carlo algorithm essentially hunts around the pa- 
rameter space seeking to reduce \ 2 to an acceptable value. A set of input parameters is 
declared to be a solution if \y^ xp — yi\ < cr, for all i. 

The relative uncertainties associated with the charged lepton masses are quite small. 
Furthermore, the charged lepton mass matrix only depends on G and H (see Eq. (jllj) ). 
These two factors make it convenient to split the search algorithm into two phases, with 
the first phase searching for Yukawa matrices G and H that yield acceptable charged lepton 
masses and the second phase searching for a Yukawa matrix F that results in acceptable 
neutrino masses and mixings. Sometimes more than one acceptable matrix F is found for a 
given pair of matrices G and H. In such cases the sets of input parameters are considered 
to be separate solutions, since in general they yield different neutrino mass matrices. 



B. Masses, mixings and phases for vr = 20 TeV and vr = 50 TeV 

In this subsection we summarize our results for neutrino masses, mixing angles and phases 
for two choices for the right-handed scale, vr = 20 TeV and vr = 50 TeV. We also include 
some comments on m^, the effective neutrino mass for neutrinoless double beta decay. In 
the following subsection we discuss some other phenomenology of the model. 

Figure^shows the neutrino masses obtained for the case v r = 20 TeV. The data were gen- 
erated using the Monte Carlo algorithm outlined in the previous subsection. Each particular 



9 For the charged leptons we adopt relative uncertainties of 5 x 10 4 , which are larger than the experimental 
uncertainties 36]. This is done for the sake of the efficiency of our Monte Carlo algorithm. 
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FIG. 1: Light (left) and heavy (right) neutrino masses for vr = 20 TeV. For each set of three 
masses, the lightest mass is indicated by a triangle, the intermediate one by a circle and the 
heaviest by a dot. 

set of "input" parameters (Yukawa couplings and Higgs VEVs) yields three light neutrinos 
and three heavy neutrinos. The plot on the left shows the light neutrino masses and indi- 
cates that the model tends to favour non-degenerate (as opposed to quasi-degenerate) light 
neutrinos. The plot on the right contains the results for the heavy neutrinos. Approximate 
expressions for the heavy neutrino masses are given in Appendix ^ As noted there, the two 
lightest right-handed neutrinos, mf and , both have masses of order €Vr, while mf has 
a mass of order vr. This scaling is evident in Fig. ^ Even though the right-handed scale is 
20 TeV, it is not uncommon to have mf below 5 TeV. The mass splitting between mf and 
is typically of order € 2 vr, as is shown in Appendix [T] 
The plots in Fig. El show the mixing angles for the left- and right-handed MNS matrices 
as well as the mass-squared differences for the light neutrinos. The top two rows of plots 
show explicitly that the constraints on mass-squared differences and mixing angles in Table 
are indeed satisfied by the model. The bottom row shows the right-handed mixing angles 
favoured by the model. Appendix ^ contains approximate expressions for each of the right- 
handed mixing angles, noting that 

sin 2 #f 2 ~ 0.5(1 -0(e)) , 

sin 2 ^ 3 ~ l-0(e 2 ), (45) 
sin 2 #f 3 ~ £>(e 4 ) . 

The interested reader is referred to this appendix for explicit expressions in terms of the 
relevant Yukawa couplings. The above expressions are consistent with the results indicated 
in Fig. H 

Figure El shows relations among the various phases appearing in the left- and right-handed 
MNS matrices, as well as plots of vl and mpp. The upper left plot shows the correlation 
between a and 6l (the phases associated with the bidoublet and left-handed triplet Higgs 
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FIG. 2: Frequency plots of mixing angles and mass-squared differences for vr = 20 TeV. The 
top and bottom rows show results for the left- and right-handed mixing angles, respectively. The 

the mass-squared differences for the light 



middle row shows frequency plots for Am^ 



neutrinos. The plots in the top two rows satisfy the constraints noted in Table [IJ 



boson fields, respectively). It is evident from the plot that the model (or at least the 
numerical procedure) favours a near and tt, although other values for a are not ruled out. 
The middle plot in the top row shows the values obtained for vl- As expected from Eq. (|35j). 
vl is of order 0.1 eV for vr = 20 TeV. The upper right plot in the figure shows that mpp, 
the effective neutrino mass for neutrinoless double beta decay (see Eq. ((3210, * s typically of 
order 0.001 or 0.002 eV for vr = 20 TeV in this model. 10 Such values are probably beyond 
the sensitivity of neutrinoless double beta decay experiments of the near future. To see why 
mpp is so small, consider its dependence on the left-handed Majorana phases a,\ and a 2 and 
the "Dirac" phase 5l (shown in the middle pair of plots in Fig. |3J). To a good approximation 
(i.e., taking cf 3 ~ 1), Eq. (|32j) may be written as follows, 



Tripp 



m 3 s 2 13 



mic? 2 



m 2 sj 2 e i{a2 - ai - n) ) e i{ai+25L) 



(46) 



As is evident from Fig. |SJ to a good approximation, a 2 ~ ai + vr (mod 2tt) and there 
is typically a partial cancellation between the terms proportional to rrii and m 2 in the 



10 This plot may be compared with Fig. 1 in Ref. |33( (although slightly different experimental ranges were 
used for that plot). 
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FIG. 3: Plots of various phases, vl and mpp for vr = 20 TeV. The effective neutrino mass mpp is 
defined in Eq. (1321 . 

above expression. Also, the m 3 term is suppressed by sin 2 # 13 . There could in principle be 
interference between the terms, but the overall smallness of would make it difficult to use 
mpp as a probe of the phases involved. The bottom row of plots in Fig. |3] shows the phases 
associated with the right-handed MNS matrix (see Eq. (|31|1). The right-handed Majorana 
phases af and satisfy the approximate relation ~ af + n (mod 2n), a result that has 
been derived analytically (see Appendix [T]) . Scatter plots for the other right-handed phases 
are also shown. 

A similar analysis has been performed for v r = 50 TeV and the results are qualitatively 
similar to those for vr = 20 TeV. The light and heavy neutrino masses obtained for that 
case are shown in Fig. 0] As might be expected, the heavy neutrino masses are larger than 
their counterparts for v r = 20 TeV and the light neutrino masses are somewhat smaller, due 
to the two seesaw mechanisms at work (see the discussion around Eq. (02))). Plots of the 
mixing angles and phases for vr = 50 TeV are similar to those in Figs. |U and El and are not 
shown. The values obtained for vl and m^p are generally somewhat smaller for v r = 50 TeV 
than those that were found for vr = 20 TeV. 
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FIG. 4: Same as in Fig. ^ but for vr = 50 TeV. Note that the axes are different here. 



C. Some phenomenological implications 

In the previous subsection we discussed neutrinoless double beta decay in the context of 
this model. Let us now consider some other phenomenological consequences of this model, 
again taking vr to be in the range 20 to 50 TeV. Incidentally, while we have not been much 
concerned in this work with the quark sector of the theory, we should note that many authors 
have studied hadronic consequences for a right-handed scale in the several- TeV range, such 
as effects on B-B and K-K mixing (see, for example, Refs. || 0, and references 

therein.) We shall not consider such effects further here, but shall be mainly concerned with 
leptonic phenomenology. 

First let us consider the effect that a "moderate" right-handed scale has on leptogenesis. 
Leptogenesis provides a mechanism for generating the baryon asymmetry of the universe 
through CP asymmetries involving leptons [ill Il2l. Il3j|. Within the LRM, the asymmetries 
can occur in the decays of heavy (right-handed) neutrinos to charged leptons and Higgs 
bosons, as well as in the decays of left-handed Higgs triplets to pairs of charged leptons (see 
Refs. QMS). 

The asymmetries arise through the interference of the tree-level diagrams 
with one-loop self-energy and vertex correction diagrams. The asymmetries for the decay 
of the lightest right-handed neutrino may be separated into Type I and II contributions as 
follows [3 EE E3 (the reason for the "Type I" designation for the first expression will be 
more apparent in a moment), 
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,(47) 
(48) 
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where Mlr = MlrVr (with Vr being the matrix that diagonalizes Mrr - see Eq. (|2*Uj) ). 

= (m^/mf) 2 and y = (m^ L /mf) 2 , with mi t being the left-handed Higgs triplet 



Also, 



mass (see Ref. |10|). It is instructive to consider the limits Xj 3> 1 and y 3> 1, in which case 
the asymmetries become 
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(49) 
(50) 



illustrating a nice symmetry between the two expressions 15]. (Recall that Ml and M„ are 
the Type I and II contributions to the light neutrino mass matrix - see Eq. (J2J.) We may 
use the above expressions to estimate the CP asymmetries within the context of this model. 
Assuming phases of order unity, ~ 10-25 TeV and M 1 / 1 ~ 0.05 eV (from Figs. □ and 



HJ), one obtains the estimates e l VR , e VR ~ 10 12 . Unfortunately, such values are far too 
small to account for the observed bar yon asymmetry of the universe - one typically requires 
e UR to be of order 10~ 6 or 10 -7 (lil l33|. We have computed the asymmetries numerically 
for the vr = 20 and 50 TeV data sets using the original expressions in Eqs. (J37j) and 
and setting m^ L = Vr for simplicity. We find numerically that 



10 



-13 



with values 



sometimes of order 10 . The asymmetry 



'"Hi 



is typically of order 10 12 to 10 11 , but 



is sometimes enhanced by one or more orders of magnitude. 11 The root cause of the tiny 
asymmetries is the fact that the right-handed scale is so low - Eqs. (J4U|) and (|5U|) are both 
proportional to the lightest right-handed neutrino mass, which is in turn proportional to the 
right-handed scale. If one were to consider a much higher right-handed scale (while keeping 
M v fixed at its physical value of approximately 0.05 eV), one could obtain asymmetries that 
are of the correct order of magnitude for leptogenesis. 

While leptogenesis would require a much higher right-handed scale than we are consider- 
ing in this work, a low or moderate right-handed scale has the phenomenological advantage 
that departures from the SM could be observable at upcoming experiments. One striking 
experimental signature of the LRM would be the production of like-sign leptons due to the 
decay of doubly-charged Higgs bosons, S LR . Several authors have investigated the possibil- 
ity of producing doubly-charged Higgs bosons at upcoming collider experiments such as the 



LHC [iol . 41 . 42j ] or a linear collider 43[ 44]|. We shall not consider direct Higgs production 
further here, except to note that a lower right-handed scale is obviously desirable if one 
hopes to produce on-shell, doubly-charged Higgs bosons. The doubly charged Higgs bosons 
in the LRM also generically contain lepton flavour violating (LFV) couplings, which are 
related to the complex symmetric matrix F in Eq. (JBJ). These couplings lead to decays such 
as /i — ► 3e [45], r — > 3fi and r — > nee. Since the decays occur at tree-level in the LRM, 
they may be used to place indirect limits on various combinations of the LFV couplings 
and the doubly- charged scalar masses. We will consider some of the limits for these decays, 
as well as branching ratio predictions for various combinations of LFV couplings and Higgs 



11 The approximation Xj 3> 1 is not a very good one for this model since X2 — 1 + 0(e), as may be inferred 
from Eqs. I|A.17|) and l|A.20|) in the Appendix. The factor 1/(1 — Xj) in Eq. I|47|l diverges for X2 — > 1. 
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TABLE II: Doubly-charged Higgs masses and corresponding branching ratios for various LFV 
decays. As described in the text, we assume that m/± L = m& R = itia, that Fl and Fr scale in 
approximately the same way as F (see Eq. (|36|)) and that the "order unity coefficients" in these 
matrices each take on the value "3." The entries in the third column give Higgs boson masses that 
would yield branching ratios at the current experimental limits. The branching ratios for the rare 
r decays are taken from Ref. and that for — * e + e~e~ is taken from Ref. The sixth 

column gives the branching ratios that would be obtained for the Higgs boson masses indicated in 
the fifth column. 
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boson masses. Other bounds on the elements of F can also be obtained by considering 



muonium-antiuonium conversion |43l. I45| and fi — > pH 47, 

Since the 2-2 and 2-3 elements of F are generically rather large in our model (see Eq. (|36|). 
let us first consider the decay r~ — > Neglecting the muon masses compared to m T , 

we obtain the following expression for the partial width, 
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\i /i jj, 



2 (32)(192tt 3 ) 
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(51) 



where the factor of 1/2 accounts for the presence of two identical particles in the final state 
and 



L,R 



v£%fvI >r 



(52) 



The matrices V[ R in the above expression are those that were used to diagonalize the charged 
lepton mass matrix in Eq. (|12|1. V[ and Vr both typically involve small mixing angles (see 
the discussion in Appendix QJ . A cross-term involving Fl and Fr has been dropped from 
Eq. (JoTj) because its leading behaviour is proportional to m 3 . The masses of the doubly- 
charged Higgs bosons were calculated approximately in Ref. ^(j for this model. For our 
purposes it is sufficient to keep the leading terms, 



m, 



2 P3 



Pi VR ' 



2p 2 v 



R > 



(53) 
(54) 



where the constants pi, p 2 and are dimensionless coefficients in the Higgs potential that 
are generically of order unity in the horizontal symmetry model. 

The first row of Table |H] gives a few numerical estimates for the decay t~ — > fi + fi~p~ 
assuming, for the sake of simplicity, that m& L = m/\ R = rriA We also assume that Fl and 
Fr both scale in the same way as F (see Eq. ([36))) and that the "order unity" coefficients in 
the elements of Fl and Fr all have a magnitude of "3." Under these assumptions, the left- 
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and right-handed contributions to Eq. (foTj) are equal. Results are also included for several 
other LFV decay modes. The expressions for those decays are similar to Eq. ([51)1 . with 
appropriate changes in the matrix elements of Fl and -Fr 12 and the deletion of the factor of 
"1/2" in cases in which the final state does not contain identical particles. The entries in 
the fourth column of the table give upper limits on the experimental branching ratios at the 
90% c.l. The third column lists the Higgs boson masses that would yield branching ratios 
right at the experimental limits. The decay /i — > 3e currently gives the furthest reach in 
terms of the doubly-charged Higgs boson mass, even though Fn is generically the smallest 
element of F in this model. The fifth and sixth columns list some representative values for 
slightly larger Higgs boson masses and what the corresponding branching ratios would be. 
For the rare fi decay we have assumed an order of magnitude improvement in the sensitivity. 
For the rare r decays we have assumed a branching ratio of 2 x 10 -9 , consistent with the 
"several times 10~ 9 " sensitivity possible at a Super B factory [50]. A Super B Factory could 
probe doubly-charged Higgs boson masses at the 20 TeV level. 

IV. DISCUSSION AND CONCLUSIONS 

In this work we have studied the lepton sector of a Left-Right Model with a low right- 
handed symmetry breaking scale. Such a model can be made phenomenologically viable if 
the underlying SU{2) L x SU{2) R x £7(1) gauge symmetry of the theory is supplemented by 
a (7(1) horizontal symmetry to suppress relevant Yukawa couplings. We have analyzed the 
parameter space of this model numerically with the help of a Monte-Carlo approach. The 
model is able to reproduce the main features of the lepton mass spectrum and is able to 
accomodate experimental constraints on the mixing angles and mass-squared differences of 
light neutrinos. 

We have also discussed other phenomenological applications of this model, such as lepton- 
flavor- violating transitions, which occur due to the presence of doubly-charged Higgs bosons 
in the theory. These transitions produce the most striking experimental signatures of the 
model. We have considered LFV decays of the type r ± — > // , which constrain both 
the masses and couplings of these Higgs bosons. While the currently-available experimental 
bounds on such decays probe the mass ranges of a few TeV for the doubly charged Higgs 
bosons, a two-order of magnitude improvement in the experimental bound for the branching 
ratios of — > fj^^/ji^ and r ± — > e T /i ± / u ± would probe the relevant mass scale of tens 
of TeV. We have noted that a LRM with such a low right-handed symmetry breaking 
scale does not accomodate the required CP- violating asymmetries needed for generating the 
baryon asymmetry of the universe via leptogenesis. Leptogenesis would generally require a 
much higher right-handed scale. We have also calculated the effective mass for neutrinoless 
double beta decay for this model. The resulting values are probably beyond the sensitivity 
of experiments of the near future. 



12 Interestingly, the Feynman rule for the vertex liL,RijL,R^\\ 1S proportional to FL : Rij = Fl,rh (no factor 
of "i"), both for the case i — j and the case i ^ j. In both cases there are two distinct contributions to 
the amplitude; these contributions are equal and cancel the factor of i. 
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APPENDIX: DERIVATION OF SOME APPROXIMATE RESULTS 

In the first part of the Appendix we derive some approximate results for right-handed 
neutrinos. Following that we comment on the possibility of obtaining large mixing angles in 
this model. 

1. Approximate relations for right-handed neutrinos 

As is clear from Fig. |2l the mixing angles for the right-handed MNS matrix are such that 
sin 2 #^ ~ 0.5, sin 2 9^ 3 ~ 1 and sin 2 ^ ~ 0. In this subsection of the Appendix we outline 
an approximate calculation of the three mixing angles, showing that, indeed, 

sin 2 #f 2 ~ 0.5(1 -0(e)) , 

sin 2 9* * l-0(e 2 ), (A.l) 
sin 2 #f 3 ~ 0(e 4 ) . 

More precise expressions for the right-handed mixing angles may be found below. We also 
determine approximate analytical expressions for the heavy neutrino masses. 

To determine the right-handed MNS matrix, one must first diagonalize the charged lepton 
mass matrix, Mg, and the mass matrix for heavy (right-handed) neutrinos, Mrr. Diagonal- 
ization of these two matrices yields the unitary matrices V[, and V£ (see Eqs. (J12)) and 
PU])). The latter two of these are used to form the right-handed MNS matrix, as seen in 
Eq. 

u% NS = vgvz . 

The "tilde" indicates that the MNS matrix still needs to be rephased to bring it into the 
usual form - see Eqs. and (|24|). Our approximate calculation makes use of the fact 
that we have two small quantities with which to perform an expansion. The first, e, is the 
parameter that breaks the horizontal symmetry in our model. Although this parameter is 
not particularly small (we take e = 0.3 throughout this paper), it still does allow for some 
progress. The second small parameter in our calculation is the ratio of the bidoublet Higgs 
VEVs, assumed to be k 2 /ki = 3/181 ~ e 4 . 

Consider first the charged lepton mass matrix, given in Eq. (jllj) . Taking into account 
the small value for the ratio A^/fci and the scaling of the Yukawa matrices in (|36|h it is clear 
that Me is dominated by the term proportional to hi, i.e., 

M e ~ —j=Hk x . (A.2) 
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(The term proportional to k 2 is suppressed by an overall factor of approximately e 8 .) Since 
H is Hermitian and k\ is real, Mi is also Hermitian in this approximation. Thus, to a good 
approximation, Mi may be diagonalized by a single unitary matrix (c.f. the general biunitary 
transformation shown in Eq. ()12|)). To the extent that Mi is Hermitian, the resulting 
diagonalized mass matrix has real eigenvalues, although some of them would possibly be 
negative. A diagonal sign matrix (±1 along the diagonal) can then be used to correct the 
signs on the masses. Thus, to a good approximation, V R and V[ in Eq. (JT2J) are related, 

V[ ~ V^A s r , (A.3) 

where A s e lg,n is a sign matrix used to make the charged lepton masses positive, and we have 



Mf ag ~ Af^Vg -^Hki V£ . (A.4) 




It is convenient to parameterize V R as a product of three 2x2 unitary rotations, as in 
Eq. (j3Ti|) . For our purposes we define 

Vi = At 6'), (A.5) 

where At is a diagonal phase matrix. The diagonalization proceeds in the following order: 
(a) diagonal phase rotation, (b) orthogonal 2-3 rotation, (c) unitary 1-3 rotation (including 
the phase 5 e ), (d) orthogonal 1-2 rotation. The form of the matrix H, 



fe w e 8 e 7 



e 8 e 6 e 5 , (A.6) 



V e 7 e 5 e 4 



H 



allows us to treat all of the 9\a in Eq. (jA.5|) as small quantities. The approximate diago- 
nalization yields 9[ 2 ~ 0(e 2 ), 9^ ~ 0(e) and 9[ 3 ~ 0(e 3 ) (we do not give the analytical 
expressions here). 

A similar procedure may be followed to diagonalize the matrix Mrr = Fvr/\^2. Since 
Mrr is complex-symmetrix (not Hermitian), the diagonalization is performed using a unitary 
matrix and its transpose (rather than a unitary matrix and its Hermitian conjugate), 

rdiag 




MT S = Vr ^Fv r V» r . (A.7) 



We parameterize V R as 



V% = A v z/°)(%, 8 V )B V , (A.8) 



where A v and B u are both diagonal phase matrices. The form of the complex-symmetric 
Yukawa matrix F, 

/e 4 

(A.9) 

yields some clues as to how to proceed with the diagonalization. First of all, we wish to order 
the eigenvalues in ascending order (by magnitude), so the 2-2 element needs to be moved 
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to the 3-3 location. This is accomplished using the 2-3 rotation that occurs immediately 
following the initial phase rotation. Since 6^ is evidently close to 7r/2, we take cos 9% 3 to be a 
small quantity for the purpose of our approximate calculation. Q v 1?> is similarly regarded as a 
small quantity. Performing these first three rotations (the diagonal phase rotation associated 
with A u , the 2-3 rotation and the 1-3 rotation) yields the following partly diagonalized mass 
matrix, 



M, 



RR 1 23,13 




(A.10) 



Clearly the remaining 1-2 rotation must be "large." Neglecting the 1-1 element in the above 
expression, it is straightforward to determine an approximate expression for sin Q v 12 in terms 
of the 1-2 and 2-2 elements and to verify that Q\ 2 ~ it/ 4. 

Combining the approximate expressions obtained for Vr and , we may finally determine 
an approximate expression for the right-handed MNS matrix. Expanding the expression in 
terms of e, we obtain the following approximate relation for sin#f|, 



sin 6^ ~ 



sin 6 V 12 



4|a| 



y/2\t 



W[\b\ + J4\a\> + \b\ 



1/2 



(A.ll) 



where 



a = F 13 = 0(e) 



(A.12) 



and 



(A.13) 
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Then, 



sin0f 2 =^= (1-0(6)) 



For the other two mixing angles we obtain 



cos 6^ 



sin 6 



13 



-^23 _|_ -^23 
-^22 -^33 



0(e) 



± 12 

-^22 



-^1*3-^23 
-^2 2 2 



H12H33 — H^H. 



13-^23 ) 



(-^22-^33 
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0(e 2 



(A.14) 

(A.15) 
(A.16) 



(Recall that F 2 2, Hi 2 , Hi 3 , H 33 and H 2 2 are all taken to be real in this work, as discussed 
in Sec. IIIBI ) Note that the 2-3 and 1-3 mixing angles receive comparable contributions 
from the diagonalizations of Mi and Mrr, while the 1-2 mixing angle is determined almost 
exclusively by the diagonalization of Mrr. The above approximate expressions for the three 
right-handed mixing angles agree relatively well with the results obtained by performing the 
diagonalizations numerically, except for cases in which the approximations break down (such 
as when a denominator in one of the expressions involved is accidentally close to zero). 
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Having performed an approximate diagonalization of M RR to obtain the mixing angles, 
it is also straightforward to obtain approximate expressions for the three heavy neutrino 
masses. We find 



mf ~ 



m 



vr 



4|a| 2 + |6| 2 - \b\ 

^\a\ 2 + W + \b\ 
22 = 0{v R ) , 



0(ev R ) 
0(ev R ) 



(A.17) 
(A.18) 
(A.19) 



which leads to the following approximate relation between mf and mf in this model, 

O(eV) • (A.20) 



m 



R 



v r 
72 



F 2 

TP 23 
^33 - -TT- 

-T22 



We may also determine the largest value we might expect for mf in our numerical work 
(given that "order unity" coefficients are required to have a magnitude between zero and 3). 
Setting 6 = and \a\ = 3e we obtain 



m 



l,max 



3ev 



R 



V2 



(A.21) 



For vr = 20 TeV, mf max ~ 13 TeV and for v R = 50 TeV, mf max ~ 32 TeV, consistent 
with Figs. ^ and 01 respectively. The largest values for mf typically occur when F 2 2 ~ 3, 
yielding m£ max ~ 42 TeV for v R = 20 TeV and m^ max ~ 110 TeV for v R = 50 TeV. The 
approximate expressions for the masses of the three heavy neutrinos agree relatively well 
with the values obtained numerically (except for cases in which the approximations break 
down, as described above). 

The approximate diagonalization procedure also allows us to derive an approximate rela- 
tion between ctf and af , two of the right-handed Majorana phases. To a good approximation 
we obtain af ~ ctf + 7r (mod 27r), a relation that is evident in Fig. El 



2. A few comments on obtaining large mixing angles 

The authors of Ref. [Hf noted that it would be difficult or impossible to obtain large 
1-2 mixing in this model. Looking at Eq. (J42|) . it would appear that the 1-2 and 2-3 angles 
should indeed be small^genericalfy. Yet our Monte Carlo algorithm does find order-unity 
sets of coefficients for F, G and H that yield neutrino masses and mixings consistent with 
experiment (i.e., with large 1-2 and 2-3 mixing angles). To understand how this can happen, 
note that the diagonalization of the neutrino mass matrix may be understood to proceed in 
several steps, beginning with a 2-3 rotation, which is followed by a 1-3 rotation and finally 
by a 1-2 rotation. 13 Experimentally, the 2-3 rotation is large (of order 7r/4), the 1-3 rotation 
is small and the 1-2 rotation is large (of order 7r/6 or 7r/5). To simplify our discussion, let us 



For simplicity we will ignore the diagonalization of the charged lepton mass matrix in our considerations 
here and consider only the neutrinos' contributions to U^ INS (see Appendix IT|l . 
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consider a case study of a neutrino mass matrix that was obtained for the case vr = 20 TeV. 
The magnitudes of the elements of the mass matrix had the following values after each stage 
of the diagonalization: 

/ 0.0015 0.010 0.011 \ 

0.010 0.0070 0.0021 eV 
\ 0.011 0.0021 0.052 J 

( 0.0038 0.011 \ 

0.011 0.0070 eV 
\ 0.054 J 

( 0.010 \ 

0.014 eV. 
\ 0.054 J 

The mixing angles for the diagonalization of this mass matrix were "typical;" i.e., the 2-3 
rotation angle was approximately it/ (3.6), the 1-3 rotation was small and the 1-2 rotation 
was approximately tt/ (5.6). The first thing that is clear is that the original mass matrix does 
not have the hierarchy described in Eq. (|42|). In particular, the 2-3 block is approximately 
degenerate and is suppressed relative to the "generic" expression. That the 2-3 block is 
quasi-degenerate is perhaps not a surprise, given that the 2-3 rotation angle needs to be 
close to 7r/4. The elements of the 1-2 block of the original mass matrix have approximately 
the expected magnitudes, with the ratio \M vl2 /M u2 2\ being relatively small. Following the 
(large) 2-3 rotation and the (small) 1-3 rotation, the 1-2 block is in a form suitable for a 
relatively large 1-2 rotation. The neutrino mass matrix in this case study is fairly typical 
in that the neutrino mass matrices that yield experimentally viable mixing angles tend to 
have suppressed (and quasi-degenerate) 2-3 blocks compared to the "generic" expectation in 
Eq. (J32J)- Some of the suppression is due to a suppression of F _1 compared to the expression 
in Eq. (JUJ). 



/ 0.0015 0.0061 0.014 
0.0061 0.033 0.025 
\ 0.014 0.025 0.020 



eV 
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